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In a partially observed quantum or classical system the information that we cannot access results 
in our description of the system becoming mixed, even if we have perfect initial knowledge. That is, 
if the system is quantum the conditional state will be given by a state matrix p r {i), and if classical, 
the conditional state will be given by a probability distribution P r (x,i), where r is the result of the 
measurement. Thus to determine the evolution of this conditional state, under continuous-in-time 
monitoring, requires a numerically expensive calculation. In this paper we demonstrate a numerical 
technique based on linear measurement theory that allows us to determine the conditional state 
using only pure states. That is, our technique reduces the problem size by a factor of JV, the 
number of basis states for the system. Furthermore we show that our method can be applied to 
joint classical and quantum systems such as arise in modeling realistic (finite bandwidth, noisy) 
measurement. 

PACS numbers: 03.65.Yz, 42.50.Lc, 03.65.Ta 



I. INTRODUCTION 

To obtain information about a system, a measurement 
has to be made. Based on the results of this measurement 
we assign to the system our state of knowledge. For a 
classical system this state takes the form of a probability 
distribution P(x' ' , í), while for a quantum system we have 
a state matrix p(t). [30] In this paper we are concerned 
with efficient simulation tcchniques for partly observed 
systems; that is, systems for which the observer cannot 
obtain enough information to assign the system a pure 
state, P(x',t) =S[x' - x(t)} or p{t) = \ip(t))(ip(t)\. 

The chief motivation for wishing to know the condi- 
tional state of a system is for the purpose of feedback 
control [1, 2, 3, 4]. That is because for cost functions 
that are additive in time, the optimal basis for control- 
ling the system is the observer's state of knowledge about 
the system. Even if such a control strategy is too difhcult 
to implement in practice, it plays the important role of 
bounding the performance of any strategy, which helps 
in seeking the best practical strategy. 

It is well known that the quantum state of an 
open quantum system, given continuous-in-time mea- 
surements of the bath, follows a stochastic trajectory 
through time [5] . In the quantum òptics community this 
is referred to as a quantum trajectory [6, 7, 8, 9, 10, 
11, 12, 13, 14, 15]. The form of this trajectory can be ei- 
ther jump-likc in nature or diffusive depending on how we 
choosc to measure the system; that is, the arrangement of 
the measuring apparatus. In this paper we review quan- 
tum trajectory theory for partially observed systems by 
presenting a simple model: A three level atom that emits 
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into two separate environments, only one of which is ac- 
cessible to our detectors. Such partially observed systems 
cannot be described by a stochastic Schròdinger equation 
(SSE) [6, 7, 8], but rather requires a more general form 
of a quantum trajectory that has been called a stochastic 
master equation (SME) [9]. This is an instance of the fact 
that the most general form of quantum measurement the- 
ory requires the full Kraus representation of operations 
[16, 17], rather t han just measurement operators [17]. 

It is also well known that if we have a classical system 
and we makc mcasurements on it with a measurement ap- 
paratus that has associated with it a Gaussian noise, then 
the evolution of this classical state in the continuous-in- 
time limit obeys a Kushner-Stratonovich equation (KSE) 
[18]. To review these dynamics for partially observed 
classical systems we present the KSE for a system that 
experiences an 'internal' unobscrvablc white noise pro- 
cess. That is, the evolution in the absence of the mea- 
surements is given by a Fokker-plank equation [19]. This 
is the classical analoguc to the quantum master equation. 

The new work in this paper is a simple numerical tech- 
nique that allows us to reduce the numerical resources 
required to calculate the continuous-in-time trajectòries. 
This method relies on the implcmcntation of linear or 
'ostensible' [12] measurement theory, classical [18] and 
quantum [10, 11, 12, 13]. For the classical case our 
method reduces the problem from solving the KSE for 
the probability distribution to simulating the ensemble 
average of two coupled stochastic diffcrcntial equations 
(SDE). For the quantum case our method reduces the 
problem from solving a conditional SME to simulating 
the ensemble average of a SSE plus a c-number SDE. 
Thus in both the classical and the quantum case, our 
method reduces the size of the problem by a factor of 
N, the number of basis states required to represent the 
system. 



2 



Recently Brun and Goan [20] havc uscd a similar 
idea to investigate a partially observed quantum system. 
Howcver, since they did not use measurcmcnt thcory with 
ostensible probabilities, their claim that they can gen- 
erate a typical trajectory conditioned on some partial 
record R is not vàlid. This is demonstrated in detail in 
A. (In their method, the record R can only be generated 
randomly, and can be found only by doing the ensemble 
average over the fictitious noise, but that is not the issue 
of concern here.) 

Finally, we combine these theories to consider the fol- 
lowing case: a quantum system is monitorcd continu- 
ously in time by a classical system but we only have 
access to the results of non-ideal measurements per- 
formed on the classical system. Note that such joint 
systems have recently been studied by Warszawski et al 
[21, 22, 23] and Oxtoby et al [24]. Warszawski et al con- 
sidered continuous-in-time monitoring of a quantum op- 
tical system with realistic photodections while Oxtoby et 
al considered continuous-in-time monitoring of a quan- 
tum sòlid- state system with a quantum point contact. 
We show that our ostensible numcrical tcchnique can be 
applicd to these types of systems, greatly simplifying the 
simulations. 

The format of this paper is as follows. In Secs. II and 
III we review quantum and classical measurement the- 
ory rcspectively. This is essential as it allows us to de- 
fine both the notations and the physical insight that will 
be used throughout this paper. In Secs. IV, V, and VI 
we investigate the above mentioned quantum, classical 
and joint systems respectively, and present our ostensi- 
ble numerical technique for each specific case. Finally in 
Sec. VII we conclude with a discussion. 



II. QUANTUM MEASUREMENT THEORY 
(QMT) 

A. General theory 

In quantum mechanics the most general way we can 
represent the state of the system is via a state matrix 
p(t). This is a positive semi-definite operator that acts 
in the system Hilbert space Tí s - In this paper we take the 
view that this represents our state of knowledge of the 
system. Taking this view allows us to simply intèrpret 
the "collapse of the wavefunction" , upon measurement, 
as an update in the observer's knowledge of the system 
[25, 26]. If we now assume that we have a measurement 
apparatus that allows us to measure observable R of the 
system, then the conditional state p r (t') of the system 
given result r is determined by [16] 



pJt') 



Òr(f,t)p(t) 

P(r,t>) ' 



(2.1) 



time. Here Ó r (t',t) is known as the operation of the 
measurement and is a completely positive superoperator 
and for efficient measurements can be defined by 

Ò r {t',t)p(t) = J[M r (T)]p(t) = M r {T)p{t)Ml(T), 

(2.2) 

where M r (T) is called a measurement operator. The 
probability of getting result r is given by 

P(r,t') = Tr[Ò r (t',t)p(t)} = Tr[F r (T)p(t)}, (2.3) 

where the set {F r (T) = M^(T)M r (T)} is the positive 
operator measure (POM) for observable R. By complete- 
ness, the sum of all the POM clements satisfics 



(2.4) 



So far we have only considered efficient, or purity- 
preserving measurements. That is if p(t) was initially 
|^i(í))(^(í)| then the state after the measurement would 
also be of this form. In a more general theory we must 
dispense with the measurement operator M r (T) and de- 
fine the Kraus operator K r j(T) [16]. This has the effect 
of changing the defmition of the operation of the mea- 
surement O r (t',t) [Eq. (2.2)] to 

Ò r {t?,t)=Y,J[KrAT% (2-5) 
/ 

and the POM elements for this measurement are now 
given by 



F r {T)=Y,Klj{T)K rJ (T). 
f 



(2.6) 



where P(r,t') is the probability of getting result r at 
time t' — t + T, where T is the measurement duration 



Note F r (T) still satisfies the completeness condition 
[Eq. (2.4)]. We can think of / as labelling results of 
fictitious measurement. 

If one is only interested in the average evolution of the 
system, this can be found via 

p(t')=Y / P(r)Pr(t') = Ò(t',t)p(t), (2.7) 



where 0(t',t) = ^2 r Ò r (t',t) is the non-selective opera- 
tion. 



B. Quantum trajectory theory 

Quantum trajectory theory is simply quantum mea- 
surement theory applicd to a continuous in-time moni- 
torcd system [6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. In con- 
tinuous monitoring, repeated measurements of duration 
T = dt are performed on the system. This results in the 
state being conditioned on a record R, which is a string 
containing the results rj, of each measurement from time 
to í but not including time 0. Here the subscript k 
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refers to a measurement complcted at time tk = kdt. 
From the record R, the conditioned state at time t can 
be writtcn as 



(2.8) 



P(R)' 

where pn(t) is an unnormalized state defined by 

|5r(í) = ò r „ (í fe , í fe _i) . . . Ó r2 (t a , ti)O ri (tx, 0)p(0). (2.9) 

The probability of observing the record R is 

P(R) = Tr[p R (í)]. (2.10) 

If we now assume that the coupling between the ap- 
paratus (bath) and the system is Markovian then the 
average state 

P (t) = ó(t k , tk.,) . . . 6{t 2 ,h)6(t x , o)p(o) (2.ii) 

is equivalent to the reduced state 

/9 re d(t)=Tr bath [|*(í))(*(í)|], (2.12) 
which itself obeys the Master equation [27] 

p(t) = Cp{t) = -i[H, P (t)} + £ lA^Mt)- ( 2 -!3) 



Here T>[À] is the superoperator defined by 
t>[À]p = Àpti - À^Àp/2 - pÀ ] À/2, 



(2.14) 



and represents dissipation of information about the sys- 
tem into the baths. 



C. Fictitious quantum trajectòries: the ostensible 
numerical technique 

If the system is only partly observed (/ in Eq. (2.5) 
represents the unobservablc processes) this state will be 
mixed. This is not a problem for simple systems but for 
a large system a numerical simulation for /?r(í) would 
be impractical. This brings us to the goal of this section 
which is to demonstrate that Pr(í) can be numerically 
simulated by using SSEs, requiring lcss space to store on 
a computer. 

To do this we assume that a fictitious measurement 
with record F is actually made on the unobservable pro- 
cess. Then we can expand /Or(í) to 

j0 R(í) = ^PR,F(í)P(P|R), (2.15) 



where 



Prf(í) = |^rf(í))(^rf(í)|. (2.16) 



Here I^R^t)) is a normalised state conditioned on both 
F and R. In quantum trajectory theory this is defined 



I^r,f(í)) 



(2-17) 



where 



|^r,f(í)) - M rkJk {dt)....M ríJl {dt)\^)). (2.18) 

Here Tk and fk are the results of the measurement oper- 
ator 



M rkJk (dt) = (r k \(fk\Ú(t k , í*-i)|0>|0>, 



(2.19) 



where |0)|0) is the initial bath state. This indicates that 
given that we have a real record R we can calculatc pr(í) 
from averaging over an ensemble of pure states IV'r.f^))- 
But as shown in A the fact that future real results are not 
necessarily independent from the current fictitious results 
mcans that we cannot generate single trajectòries with- 
out knowing the full solution. However by using quantum 
measurement theory with ostensible distributions we can 
get around this problem. 

Under ostensible quantum trajectory theory [11, 12, 
13] we can define a state, |'0r.f(í)) as, 



fe,F(í)) = 



VA(F.R) 



(2.20) 



where A(F,R) is an ostensible probability distribution. 
This is simply a guessed distribution that only has the 
requirement that it be a probability distribution and be 
non-zero when P(F, R) is non-zero. Note this state is 
no longer normalized to one and this is why we signify it 
with the bar. The true probability can be related to the 
ostensible probability by 



P(R,F) = (^ r , f (í)|Vr,f(í))A(F,R), 



(2.21) 



which is a generalized Girsanov transformation [5, 11, 12, 
13, 28]. 

Going back to Eq. (2.15) and using the above equations 
we can write pn(t) as 



where 



E F I^R,F(*))(^R.F(t)|A(F,R) 
P(R) 



(2.22) 



P(R) = 5J(^R,f(í)|W(*))A(F, R). (2.23) 

F 



Note that the sum containing A(F, R) in the above equa- 
tions simply represents the ensemble average over all pos- 
sible fictitious records. Thus we can rewrite Eq. (2.22) 
as 



Pr(í) 



Et 



|^r,f(í))(^r,f(í)| 



(^r,f(í)|^r,f(í)) 



(2.24) 
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III. CLASSICAL MEASUREMENT THEORY 
(CMT) 

A. General theory 

In this paper when considering what we call a classical 
system, we are referring to a system that can be described 
by the probability distribution P{x, t) (i.e a vector of 
probabilitics) rather than a state matrix. That is, with 
respcct to a fixed basis x the coherences (off diagonal 
elements) are always zero. If we now measure observable 
R of the system, then after a measurement which yielded 
result r, the state of the system is given by [29] 



P r (x,t) 



P(r,t\x,t)P(x,t) 



where 



P(r,t)= / dxP(r,t\x,t)P(x,t). 



(3.1) 



(3.2) 



This is known as Baves' thcorcm. Hcre P r (x,t) = 
P(x, t\r, t) is called a conditional state and represents our 
new state of knowledge given that we observed result r. 
Here we have only considered minimally disturbing classi- 
cal measurements. That is, there is no back action acting 
on the system in the measurement process. To general- 
izc Baycs' thcorem to deal with measurements which in- 
cur back action we mathematically split the measurement 
into a two stage process. The first is the Bayesian update, 
followed by a second stage described by B r (x' , t'\x, t), the 
probability for the measurement to cause the system to 
make a transition from x at time t to x' at time i' = t + T, 
given the result r. Thus for all x' , x and r 



B r (x',t'\x,t) > 0, 
dx'B r (x',t'\x,t) = 1. 



(3.3) 
(3.4) 



Now by defining the operation 

O r (x',t'\x,t) = B r (x',t'\x,t)P(r,t\x,t) (3.5) 

the conditional system state after the measurement be- 
comes 

P^^ÍÍWMM, (3.6) 



P{r,t>) 



where 



P(r,t') = j dx' J dxO r (x',t'\x,t)P(x,t). (3.7) 
Using Eq. (3.4) this can be rewritten as 

P(r,t') = J dxF r (x,t)P(x,t), (3.8) 
where F r (x,t) = P(r,t\x,t), which by definition satisfics 
^F r (x,t) = l, (3.9) 



is the classical analoguc of the POM clement. The aver- 
age evolution of the system is given by 

P(x',t') = ^P r (x',t')P(r,t') 

dxO(x',t'\x,t)P(x,t), (3.10) 



where 0(x',t'\x,t) = ^ r O r (x' ,t'\x,t) is the non- 
sclcctive operation. 

Note that for any B r [x! , t'\x, t) that satisfies Eqs. (3.3) 
and (3.4) we can rewrite it as 

B r (x',t'\x, t) = J2 SW ~ Xr,f{t')]P(f, t'\x, t; r, t), 

(3.11) 

where x r j(t') is the new system configuration x' at time 
i' given the measurement result r and extra noise / (the 
stochastic part of the back action). The parameter / 
is analogous to the fictitious measurement results in the 
quantum case. Thus the operation for the measurement 
can be written as 

O r (x',t'\x,t) = J2à[x' - x rJ (t')]P(f,t';r,t\x,t), 
f 

/ 

This is the classical equivalent of Eq. (2.5). 

In the above we have purposcly structured QMT and 
CMT so that the theories appear to be similar and as 
a general rule we will push this point of view through- 
out the rest of this paper. However, it is important to 
point out the key diffcrcnccs bctween these theories. In 
the quantum case we can always write the measurement 

operator (or Kraus operator) as M r = Ü r \f~F\- where U r 
is a unitary operator. That is we can always intèrpret a 

measurement as a two stage process, where \/Fr is rc- 
sponsible for the wavefunction collapse and the gain in in- 
formation by the observer and U r is some extra evolution 
that entails no information gain (as the entropy of the 
system is not changed by this evolution). It simply adds 
surplus back action to the system. In the classical case 
we can also write the measurement as a two stage pro- 
cess. However, the first process by definition has no back 
action; it is simply the update in the observer's knowl- 
edge of the system. Furthermore the second stage is not 
necessary unitary evolution (and as such can change the 
entropy of the system) . Thus back action in the quantum 
and classical case are physically diffcrent processes and 
one can not separate all the back action in the quantum 
case from the observer's information gain. Mathemati- 
cally speaking, the difference arises from the fact that a 
quantum state is represented by a positivc matrix, the 
state matrix, while a classical state is represented by a 
positive vector, the vector of probabilitics. 
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B. Classical trajectory theory 

To achieve continuous-in-time measurements theory 
for a classical system we simply let the measurement 
time tend to dt and extend the numbcr of consecutivc 
measurements to t/dt. Then the state of the classical 
system given the measurement record R is 



whcrc 



P(R) ' 



(3.13) 



wherc Pr(x, t) is an unnormalized state defined by 
Pr(x,í) = / dxk-x— I dxi l dx 



xO rjt (x,í|a; fc _i,í fc _i) . . .O r2 (x 2 ,t 2 \xi,t 1 ) 
xO ri (x 1 ,t 1 \x ,0)P{x ,0). (3.14) 

The probability of observing the record R is 

P(R) = / dxP R (x,t). (3.15) 



If we now assume that the noise added by the measure- 
ment apparatus is white, and the form of the back action 
is independent of the rcsults R, then the unconditional 
state 



P(x,t) = J dxk-i-.. J dx! J dx 



xO{x,t\x k -i 1 t k _ 1 ) . . .0(x 2 ,t 2 \xi ) ti) 
xO(x 1 ,t 1 \x Q ,0)P(x ,0) (3.16) 

is the solution of the Fokker Plank Equation [19] 

d t P(x,t) = -d x [A(x,t)P(x,t)\ + ±d*[D 2 (x,t)P(x,t)], 

(3.17) 

where A(x, t) determines the amount of drift and D(x,t) 
determines the amount of diffusion. 



C. Fictitious classical trajectòries: The ostensible 
numerical technique 

The bàsic principle behind this technique is that we 
assume that the unobservable proccss, F, that generates 
the back action part of the measurement is fictitiously 
simulated. To be more specific we can define 



Pr.f(x, t) = dx 



dxi / dxo 



x Jr k J k (x,t\x k -l,t k -l) ■ ■ ■ Jr 2 J 2 (x2,t 2 \xi,ti) 
xJ riJl ( Xl , h\xo, 0)P(X , 0). (3.18) 

where J r j{x', t'\x, t) is defined implicitly in Eq. (3.12). 
From this the conditional state, Pr(x,í), is given by 



Pn(x,t) = J2 p RAx,t)P(F\K), 



(3.19) 



Pn,F{x,t) 



PR,F{x,t) 

P(R,F) 



(3.20) 



But as in the quantum case this cannot be directly cal- 
culated and as a result we must use an ostensible theory. 
We definc the ostensible state by 



Pn,F(x,t) 



P-R.F(x,t) 
A(R,F) ' 



(3.21) 



and the true probability can be related to the ostensible 
by 



P(R,F) = J dxP R , F {x,t)A(R,F), 

the classical Girsanov transformation. 

Using the above we can rewrite Eq. (3.19) as 



Pn(x,t) 



E F FR, F (^,f)A(F,R) 
P(R) 



whcrc 



P(R)=Y / I dxP RtF (x,t)A(F,R). 

T7 1 J 



(3.22) 



(3.23) 



(3.24) 



As in the quantum case we can rewrite Eq. (3.23) as 



Pr(x, t) 



Ep 


P-R,F(x,t) 




E F 


/ dxP R , F (x,t) 



(3.25) 



where Pr,f(^,í) is an unnormalized pure classical state. 
That is, it is of the form P(x,t) = pu.fS[x — xr.fÍí)], 
where pr.f is the norm of the ostensible state. To show 
this we consider a system initially in the state P(x, 0) = 
pS(x — xq) then by using Eqs. (3.21) and (3.18) with 
Jr-í,U ( x ' , íi|x, 0) defined implicitly in Eq. (3.12) we can 
rewrite P ri j 1 (x / ,ti) as 

Pr 1 J 1 (x',tl)=Pr lJl (tl)6[x'-Xf 1 , ri (t 1 )], (3.26) 

which is still of the (5-function form. Here Pn,fi(ti) is 
given by 

PriJl (h) = P(f u h; r 1} 0\x , 0)p(0)/A(n, fi), (3.27) 

and Xf liri (ti) is determined by the underlying dynamics. 
That is, we can simulate the distribution by solving the 
two coupled SDEs, £r,f(í) and pr.f(í)- 



IV. A QUANTUM SYSTEM WITH AN 
UNOBSERVED PROCESS 



To illustrate a quantum system where a complete mea- 
surement can not be performed, due to some physical 
constraint, the system in Fig. 1 was considered. This 
system is a three level atom with lowering operators 
Li = |1)(3| and L 2 = |2)(3|, and decay rates 71 and 
72 respectively. 




FIG. 1: A simple system (a three level atom) which has two 
outputs due to the to lowering operators L\ and Li. 
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FIG. 2: The solution to the master equation. The íirst sub- 
plot shows P33(í) (sòlid line), pii{£) (dashed line) and pn(í) 
(dotted line). The second and third subplot show the real 
and imaginary parts respectively of pu{t) (sòlid line), p3i(í) 
(dashed line) and pz-iif) (dotted line). The fourth subplot 
illustrates the purity of this state. This is all for the initial 
condition |^(0)) = 0.4123|1) + 0.1|2) + (0.9 + 0.1i)|3) and 
71 = 0.5 and 72 = 1. 



A. Master equation 





M0,t] 



FIG. 3: A schematic representing homodyne measurement of 
one of the outputs of the three level atom. In an ordinary 
homodyne measurement the signal is coupled to a classical 
local oscillator (LO) via a low reflective beam splitter and 
then detected using a photoreceiver. 



this arrangement is an incfficient measurement we have 
to use the operation defincd in Eq. (2.5). To determino 
the Kraus operators we need to present the underlying 
dynamics in more detail. For the interaction of this sys- 
tem with a Markovian bath (and under the rotating wave 
approximation and in the interaction frame) the total 
Hamiltonian is 



H(t) = iVtT j 6(t - t')[L 1 bl(t') - Llb r (t')]dt' 

+ih^ í 5(t - t')[L 2 b\(t') - Llb f (t')}dt' . 



(4.1) 

Here b r (t) and 6/(í) are the temporal- mode annihilation 
operators for the detected (b r ) and non-detected (ò/) 
fields (baths). Since these fields are Markovian there will 
be a commutator rclationship for the fïcld of the following 
form 



With no external driving [H = in Eq. (2.13)], the 
solution of the master equation can be determined an- 
alytically. To illustrate a non-trivial solution we cal- 
culated this solution for the initial condition \ifi(0)) = 
0.4123| 1) +0.1 12) + (0.9 + 0.H)|3) and coupling constants 
71 = 0.5 and 72 = 1. This is shown in Fig. 2. In this 
figure it is observed that as time goes on, the state be- 
comes mixed. This is seen as the purity p(t) = Tr[p 2 (í)] 
of the state decays (although not monotonically) as time 
increases. This figure also shows that the state becomcs 
a mixture of the two ground states, with the ground 
state associated with the larger coupling constant being 
weighted more heavily, even though it started with less 
weight. 



B. Conditional evolution: The quantum trajectory 

In this section we consider the trajectory pn(t) which 
occurs when output L\ is monitorcd using homodyne- 2; 
detection and output L 2 is un-monitored. A schematic 
of this measurement process is shown in Fig. 3. Because 



[b i (t),b](s)}=6(t-s)5 i>j , (4.2) 

where i, j denotes cither of the two baths. This indicates 
that the field operators arc gaussian white noise opera- 
tors. Thus they obey Itò calculus and the infinitesimal 
evolution operator is [7, 15] 

Ü(t + dt,t) = exp|^/7Ï[iid-BÍ(í) -L\dÈ r (t)] 

+^fy[L 2 dÈ\(t) - LldÈ f (t)}}, (4.3) 

where dBi satisfies the commutator relation 

[dB{t) l ,dè]{t)]=dt8 lJ . (4.4) 

Thus U (t+dt, i) is an operator acting in the Hilbert space 
H s (£> Tí r (8) Uf , where H s , H r and H / are the Hilbert 
spaces for the system, detected field and non detected 
field respectively. 

Now, given that a projective measurement is made on 
bath field b r (t) and bath field is completely unob- 

served the state of the system after this measurement 
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(time dt later) is given by Eqs. (2.1) and (2.5) with 
T = dt, and thc Kraus opcrator is 

KrAdt) = (f\f(r\rÜ(t + dt,t)\0) r \0)f. (4.5) 

Here {|?')r} is the set of orthogonal states the bath is 
projected into, while {|/)/} is any arbitrary orthogonal 
basis set. For a homodyne-x measurement of bath b r (t) 
the set {|r) r } corresponds to the eigenset of the operator 
dB r (t) + dBJ.(t) [11] and the results r are the correspond- 
ing eigcnvalucs. Note we have assumed that initially the 
baths, for all thc tcmporal-modes, are in the vacuum 
state. 

After some simple rearrangement and using (rdt) 2 = 
dt, the POM elements for this measurement are of the 
form 



F r (dt) = | 
whcrc x\ = L\ 
P(r,t 



'■[l + v^nrit + dtfdtx!}, 



í[. Thus 



dt) = |(r|0)ni + rdV7i^i;tJ, 

where (ii)t = Tr[£ip(í)]. Using the fact that 
temporal-quadrature state, 



l(r|0) r | 

we can rearrange this to 
r dí" 



cxp 



( 2/dt) 



P(r, t + dt) 



exp 



7i \xi)t\ 



2/dt 



(4.6) 



(4.7) 



is a 



(4.8) 



(4.9) 



This implies that the random variable associated with 
this distribution, r(t + dt)dt, is a gaussian random vari- 
able (GRV) of mcan ^/jï (x\)tdt and variance dt. That 
is, 



r(t + dt)dt = dW{t) + dt y /j[{x 1 ) u 



(4.10) 



whcrc dW(t) is a Wiener increment [19]. 

Using the above and Eqs. (2.1) and (2.5) the stochastic 
master cquation for this system is 

dp R (t + dt) = dt( l2 V[L 2 ] + 7iD[Li] 



(4.11) 



+dW(t)^H[L 1 }/dt)p R (t), 



where Ti.[À] is the superoperator 

H[À]p = À P + pit - Tr[Àp + pÀï] f 



(4.12) 



To illustrate an example quantum trajectory, 
Eq. (4.11) was solved for a randomly chosen record 
R and the same parameters used in Fig. 2. This is 
shown in Fig. 4. It is observed that this state evolution 
is stochastic in time and becomes mixed (but not as 
mixed as the averagc evolution). It is interesting to note 
that by performing this measurement the coherence 
P12 r(í)j which was a constant of motion for the averagc 
state, becomes comparable to thc other coherence and 
does not decay with time. 
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time, t [y, 1 



FIG. 4: The solution to pn written in matrix elements. The 
first subplot shows P33,r(í) (sòlid line), p22,n(t) (dashed line), 
Pii,r(í) (dotted line). The second and third subplot show 
the real and imaginary parts respectively of Pi2,r(í) (sòlid 
line), P3i,r(í) (dashed line) and P32,r(í) (dotted line). The 
fourth subplot illustrates the purity. We have used the same 
parameters as in Fig. 2 



C. The ostensible numerical technique 

In Sec. II C we observed that the conditional evolu- 
tion of a partly monitored system could be simulated by 
assuming that fictitious measurements are made on the 
unobservable process. For this system we assume that 
a fictitious homodyne-x measurement is made on output 
L-2- Note we could have chosen any unraveling for F. 

To determine the SSE for the ostensible state |'0r i f(í)) 
[the state which we substitute into Eq. (2.24) to deter- 
mine the actual conditional evolution] we have to derivc 
the measurement operator for the combined real and fic- 
titious measurements, as well as makc a convenient choicc 
for A(F, R). Using Eq. (2.19) and the fact that we are 
performing homodyne-x measurements the measurement 
operator is 



M r . f (dt) = (/|0)(r|0)(l + V7Ï>&Li 
-Tídtllli/2 - l2 dtL\L 2 /2\ , 



72 fdtL 2 
(4.13) 



where the bath states |/) and |r) are temporal quadrature 
states acting in Hilbert spaces Tif and TL r respectively. 
To derive this we have expanded Eq. (4.3) to first order 
in dt and used the fact that (fdt) 2 = (rdt) 2 = dt. Since 
the real distribution is Gaussian (with a variance í/dt) a 
convenient choice for A(F, R) is A(F)A(R) where A(F) = 
A(f k ) . . . A(/i) and A(R) = A(r fc ) . . . A(n) with 



A(r) 




(r-A) 2 



2/dt 



(4.14) 
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A(/) 



exp 



2/dt 



(4.15) 



Here À and p are arbitrary paramctcrs. With these os- 
tensible distributions, Eq. (4.13), and Eq. (2.20), the os- 
tensible SSE is 



c#R,fW) = dt([r-\}(^L 1 -\/2) + [f-p] 

x(VtTí 2 - M/2) - \ \lxl\L x + l2 L\L 2 

-VtTALi - ^pL 2 + A 2 /4 + M 2 /4 

x|^r,f(í)). (4.16) 

Now since we are interested in calculating Pr(í) bascd 
on an assumed known real record R, we can rewrite 
Eq. (4.16) as 



dc\ = 



dci 



dc 3 



where 7 



c 3 [V7T(r - A)dt + dtX/2] - di^dWn 
+V7Ï (r - A)díA + dt\ 2 /A + dtp 2 /4]/2(4.17) 

+V7Ï (r - A)díA + dt\ 2 /A + dtp 2 /4}/2(A.18) 
c 3 [-7dí + ^/yïdWp + VtT (r - A)díA 
+díA 2 /4 + dípt 2 /4]/2, (4.19) 

7i +72- Here we have used the identity 

|^(í))= Cl |l)+c 2 |2)+C3|3), (4.20) 



and replaced fdt with dW(t) + p,dt, where dW(t) is a 
Wiener increment. 

To illustrate the convcrgcnce of our method the ensem- 
ble averagc of the above ostensible SSE for A = /1 = 
was calculatcd for n — 10 and n — 1000. To quantify 
how closely the ensemble method reproduces pn(t) we 
used the fidelity measure, which for two different quan- 
tum states is defined as 



F (Q) (t) = ^[\/VpT(í)p2(í)a/Mí) 



(4.21) 



Note this measure ranges from to 1 with indicating 
two orthogonal states and 1 indicating the same state. 
The result of this measure for the actual /9r(í) and the 
ensemble version are shown in part A of figure 5. Here 
we see that for larger ensemble size the fidelity is closer 
to one, indicating that as we increase the ensemble size 
our ostensible method approaches the actual pn(t). 

To illustrate the effect of choosing different ostensible 
distributions we considered the casc whcn A = and 



(^R,f(Q|^2 + £ 2 |^R,f(Í)) 

(V>r,f(í)|Vr,fW) 



(4.22) 



That is, the ostensible probability for the k th fictitious 
results is the true probability we would expect based on 
the past real and fictitious results up to, but not includ- 
ing the time kdt. The motivation for this choice is that 
with [i = 0, the improbable trajectòries, ones that tend 
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FIG. 5: This figure shows the fidelity between the actual pn(t) 
and our ensemble method for ensembles sizes 10 (dotted) and 
1000 (sòlid). Part A corresponds to a linear ostensible distri- 
bution while part B refers to the non-linear ostensible distri- 
bution. The same parameters were used as in Fig. 2. 



towards being inconsistent with the full real record, will 
have norms that are very small and as such have little 
contribution to the ensemble average. By contrast, us- 
ing Eq. (4.22), the improbable trajectòries are less likely 
to be generated, so avoiding useless simulations. With 
this ostensible distribution the fidelity measure was cal- 
culated for n = 10 and n = 1000. These results are shown 
in part B of figure 5. Here we see that for the smallcr 
ensemble size the fidelity is closer to one than that ob- 
served using the first ostensible case. This indicates that 
the rate of convergence for this case is greater than the 
A = /i = case. 



V. A CLASSICAL SYSTEM WITH AN 
INTERN AL UNOBSERVED PROCESS 

In this section we consider continuous-in-timc mca- 
surements with Gaussian precision of a classical system 
driven by an unobservable noise process. This for ex- 
ample could correspond to a measurement of the voltage 
across a resistor that is driven by a noisy classical current. 



A. The average evolution 

We restrict ourselves to unconditional state evolution 
described by the Fokker Plank equation Eq. (3.17). This 
cquation has as its solution a distribution that diffuscs 
and drifts though time. Using Eq. (3.10) and only con- 
sidering one interval in time we can write 

P(x',t + dt)= dxO(x',t + dt\x,t)P(x,t), (5.1) 
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which when comparec! to Eq. (3.17) implics that RHS of 
the above equation cquals 



dx[l - dtd x >A{x,t) +dtdl,D 2 (x,t)/2]5(x' - x) 
xP(x,t). (5.2) 



By introducing an arbitrary Gaussian distribution 
P(f,t + dt) with mean m(t) and variancc 1/dt, that is 



P(f, t + dt) = J^- exp 
V Zlï 



[f-™(t)]< 
2/dt 



(5.3) 



Eq. (5.2) can be rewritten as 

dfP(f, t + dt) J dx[l - dtd x ,A(x, t) - dt[f - m{t)] 

xd x ,D f (x,t) + dtd 2 ,D 2 (x,t)/2]6(x' -x)P(x,t). 

(5.4) 

By using Itò calculus and a Taylor expansion this can be 
rewritten as 

J dxE f ^S[x' -x- dtA(x,t) - dt[f(t + dt) - m(í)] 

xD(x,t)]\p(x,t). (5.5) 
where f(t + dt)dt = m(t)dt + dW{t). Thus 

0(x',t + dt\x,t) = Ef{6[x' -x f (t + dt)]}, (5.6) 

where Xf(t + dt) is determined by the following SDE 

dx F (t) = dtA[x F (t),t] + dt[f(t + dt) - m{t)]D[x F (t),i\. 

(5.7) 

Noto here we havc writtcn the SDE for the completo 
record F. 



B. Conditional evolution: The 
Kushner-Stratonovich equation 

To derive the KSE we start by deriving O r (x' ,t'\x,t) 
and P(r, t + dt). For the case when the classical measurc- 
ment has a back action that is independent of the result 
r(t + dt), the operation for the measurement is given by 

O r (x', t'\x, t) = 0(x', t'\x, t)P(r, t\x, t), (5.8) 

where 0(x',t'\x,t) = B(x',t'\x,t). Thus to derive 
O r (x',t'\x, t) we need only P(r,t\x,t). For a measure- 
ment that has a precision limited by Gaussian white noise 
it follows that 



P(r,t\x,t) = F r (x,t) = ^=cxp{-(r~x) 2 dt/2[3}, 
V2irp 

(5.9) 

where (3 is a constant characterizing the classical mea- 
surement strcngth. 



To find P(r, t) we substitute Eq. (5.9) into Eq. (3. 
This gives 



exp[-(r - x) 2 dt/2(3]P{x, t). 



P(r, t + dt) = dx 



dt 



121^ 



(5.10) 

After some simple stochastic àlgebra and using r 2 = (3 /dt 
this can be simplified to [22] 



P{r, t + dt) 



dt 



I2j^ 



exp[-(r- (x) t ) 2 dt/2f3], (5.11) 



where for the classical system (x)t = J xP(x,t)dx. From 
Eq. (5.11) the stochastic representation of r(t + dt) is a 
Gaussian random variable with mean (x)t and variancc 
(3dt. That is, 



r(t + dt) = (x) t + y/0dW(t)/dt. 



(5.12) 



With all the above information and Eq. (3.6) the con- 
ditional state at time t' = t + dt is 



P r (x',t + dt) 



dxEf 



{ó[x' -x f {t + dt)]\{l 
x) t ]dt//3}p(x,t). (5.13) 



Here we have expanded the exponentials in Eq. (5.11) 



and Eq. (5.9) to second order in dt and used 



0/dt. 



Taylor expanding the delta function and averaging over 
the f{t + dt) [using Eq. (5.3)] for each step in time gives 
the KSE 

Pnix^ + dt) = P R {x,t) + dt[x-(x) t ][r{t + dt)-(x) t ] 
xPr(x, t)/(3 - dtd x [A(x, t)Pn(x, t)} 
+ ±dtd 2 x [D 2 (x,t)P R (x,t)] (5.14) 

and (x)t becomes J xPYi{x,t)dx. In general to solve this 
equation we need to solve for all x. For some A(x, t) and 
D(x,t) this can be a rather lengthy numerical problem. 
In the following section we will present our ostensible 
tcchnique which allows us to reformulate the problem to 
solving two couplcd SDEs, at the cost of performing an 
cnscmblc average. 



C. The ostensible numerical technique 

As shown in Sec. III C, if we consider the unobservable 
process F as actually oceurring then we can simulate the 
KSE by using Eq. (3.25), and Pr.f(^, *) is determined by 
solving two couplcd SDEs. For the case when the classical 
measurement has Gaussian precision and the back action 
only depends on the white noise process /(i), we can 
rewrite P(f, t'; r, t\x, t) in Eq. (3.27) as P(f, t')P(r, t\x, t) 
where P(f, t') is given by Eq. (5.3) and P(r, t\x, t) is given 
by Eq. (5.9). Thus x R ^ F (t) becomes x F (t) and is given 
by Eq. (5.7). To find the diffcrcntial equation for pr,p(í) 
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we need to assume a form for the ostensible distribution 
A(/,r). We use A(/,r) = A(/)A(r), where 



A(r) = \j — CX P 



(r-A) 2 



(5.15) 



and A(/) is given by Eq. (4.15). Extending Eq. (3.27) to 
continuous measurements gives 

*r,f(í) = dí[rn(í)-/i][/(í + dí)-/i]p R ,F(í) 

+dt[x{t) - A][r(í + cfó) - A]p RjF (í)//3. 

(5.16) 

Thus to dctcrminc -Pr.f(^,í) wc only need to simu- 
late Eqs. (5.7) and (5.16) with R assumed known and 
f(t + dt) given by Eq. (4.15). Pa(x,t) is then deter- 
mined by Eq. (3.23). Since the theory requires Pa,F{x,t) 
to be a delta function, one might conclude that this 
mcthod is only vàlid for initial conditions of the form 
P(a;,0)p(0) = 8(x — Xq). Infact, we arc not limited to 
this case. To consider other initial conditions we simply 
choosc the initial value xq in Eq. (5.7) from the distribu- 
tion P(x,0). 



D. A simple example 

To illustrate the classical theory we consider a Gaus- 
sian measurement of a classical system that is driven by 
an an unobscrvable white noise process with m(t) = 
and drift and diffusion functions given by 



A(x,t) 
D(x,t) 



= —kx + l, 
= b. 



(5.17) 
(5.18) 



If this is the case then Pr(x,í) has a Gaussian solution 
with a mean (xji)t and variance ^r(í) given by 

d(x R ) t = dt{v Ii (t)[r(t + dt)-(x R )t]/P-k(x R }t 

+1}, (5.19) 
dv R {t) = díh^|(í)//?-2fc l / R (í) + ò 2 ], (5.20) 

and r(t + dt) = (ïr) ( +dW(t). That is, as time increases 
the measurement has the effect of reducing the variance 
but the diffusive coefficient b causes this variance to in- 
crease. The mean, however contains both the determin- 
istic evolution and a random term due the measurement. 
To illustrate this solution we have simulated Eqs. (5.19) 
and (5.20) for the case when A(x, t) = 1 — x, D = 1 and 
/3 = 1. The results of this simulation are shown in Fig. 6 
as a sòlid line. Here we see that the mean follows some 
stochastic path conditioned on the record R, while the 
variance follows a smooth function. 

To illustrate our ostensible method we use the above 
record and solve numcrically Eqs. (5.7) and (5.16) with 
A = /i = 0. The mean and variance is then found via 



E F 


xf(í)pr,f 




Ep 


ï>r,f(í) 





(5.21) 
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FIG. 6: The mean and variance of Pn.(x,t) when (3 = 1, 
A = 1 — x and D — 1 for both Pn.(x,t) calculated exactly 
(sòlid) and via the linear method for an ensemble size of 10 
000 (dotted). 



E F 


X F {t)p R ,F 




Ep 


í>r,f(í) 





- [xa 



(5.22) 



where Ef denotes an ensemble average over all possible 
fictitious records. The numerical vàlues for the mean 
and variance are shown in Fig. 6 (dotted) for an ensem- 
ble size of 10 000. To get an indication of the numeri- 
cal error in the solution from our method, the difference 
from the exact solution is shown in Fig. 7. The dotted 
line corresponds to an ensemble of 100 and the sòlid to 
one of 10 000. Here we see that the ostensible method 
solution agrees well with the exact solution and as wc 
increase the ensemble size the difference between these 
solutions decreases. To get a better indication of how 
well our method reproduces the actual Pr(x,í), we also 
calculated the classical hdelity, which is dchned by 



F (C) (Í) = / dxy/P!(x,t) y/P2(x,t). 



(5.23) 



This was calculated under the assumption that the state 
calculated via the ostensible method was also Gaussian. 
This is illustrated in Fig. 7, where we see that for the 
larger ensemble the hdelity is very close to one, implying 
that the distributions arc almost identical. 



VI. AN UNOBSERVABLE QUANTUM SYSTEM 
DRIVING A CLASSICAL SYSTEM 



In this section we consider the following situation: a 
quantum system is monitored continuously in time by a 
classical system. This in turn is measured with Gaussian 
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FIG. 7: The first and second plot show the difference between 
the mean and variance of Pn(x,t) calculated by the linear 
method and the know result for ensemble sizes 100 (dotted) 
and 10 000 (sòlid). The third plot shows the Fidelity between 
Pn.(x,t) calculated by the linear method and the know re- 
sult for ensemble sizes 100 (dotted) and 10 000 (sòlid). The 
parameters are the same as in Fig. 6. 



precision, and these are the only results to which we have 
access. This for example occurs when the signal from the 
quantum system enters a detector with a bandwidth B, 
resulting in the state of the detector being relatcd to F 
by [22] 



x(t) 



dsBexp[-B(t - s)]f(s). 



(6.1) 



Thus in a measurement that reveals x(t) with perfect 
precession [eg F r (x) = P(r,t\x,t) = S(r — x)] wc could 
determine F (the quantum signal) by inverting the con- 
volution in Eq. (6.1). But if this measurement has Gaus- 
sian precision [Eq. (5.9)] thcn wc must treat the statc 
of the detector as a classical probability distribution and 
use a mixture of CMT and QMT to describe the condi- 
tional state of the supersystem (classical and quantum 
system). To denote the supersystem we use the nota- 
tion p[x,t), where a; refers to the classical configuration 
space and p denotes an object acting on a Hilbcrt spacc. 
This has the interpretation whereby P(x,t) = Tr[p(x,í)] 
is the (marginal) classical state and p(t) — J p(x,t)dx is 
the (reduced) quantum state. For uncorrelated quantum 
and classical states, p{x,t) = P(x,t)p(t). 



A. Conditional evolution 

We denote the state of the supersystem conditioned 
on the classical result r at time t + dt as p r (x,t + dt). 
Assuming that the quantum system is not affected by 
the classical system, this can be expanded as 

p r (X, t + dt) = ^ P r (/. 1 + dt ) P r,f (x,t + dt)pf(t + dt), 

(6.2) 



where Pf(t + dt) is the state that an observer who had 
access to all the quantum information would ascribe to 
the quantum system. That is, fit + dt) can be regarded 
as really existing (with the collapse of the wavefunction 
oceurring at this level); it is just that the real observer 
does not have access to this information. The state of 
knowledge of this real observer is different from, but con- 
sistent with, that of the hypothetical observer who has 
access to F. 

In terms of the operation of the measurement, the con- 
ditional statc can be written as 



p r (x' ', t + dt) 



p r {x' ', t + dt) 
P(r,t + dt) ' 



(6.3) 



where 



.(x',t + dt) = J dx^2j r j(x',t + dt\x,t) 



xO f (t + dt, t)p(x, t)/P(f, t + dt) (6.4) 



and 



P(r,t + dt) = / dx'Tr p r (x',t + dt) 



(6.5) 



The quantum part of the operation of measurement in 
defined by Eq. (2.5) and the classical part is defined in 
Eq. (3.12) with the replaccmcnt of P(f,t';r,t\x,t) — > 
P(f,t')P(r,t\x,t) because in this system the quantum 
signal does not depend on the classical state. 

To illustrate the above we consider the case when we 
are monitoring with Gaussian precision the classical sys- 
tem defined by Eq. (6.1) which is in turn monitoring the 
x quadrature flux coming from a classically driven two 
level atom (TLA). This is the same as the system con- 
sidered in Rcf [22] and as such we will simply list the 
important cquations. The quantum part of operation is 
given by Ó f (t + dt, t) = J[M f (dt)] where 

M f (dt) = (f\0)[l-dt(iH-^ffa + ^à/2)}. (6.6) 
The fictitious quantum signal statistic obeys 



P(f,t + dt)= / dxTr[O f (t + dt,t)p(x,t)], 



(6.7) 



which for a homodyne-a; measurement can be shown to be 
of the form displayed in Eq. (5.3) with m(t) = ^/jTr[(à+ 
à^)p(t)]. Here tr is the lowering operator for the TLA 
and 7 is the decay rate. Note here we have assumcd 
all the quantum signal is fed into the classical system, if 
we wanted to simulate some inemeieney we would simply 
use the Kraus represention, and for the case where this 
inemeieney is a constant, 77, we simply replace a in the 
above equations by -y/fja. 

As shown in Sec. V for a classical measurement with 
Gaussian precision and a back action that does not de- 
pend on the results of the measurement, the classical part 
of the operation is 



J rJ {x' ,t+dt\x,t) 



5[x'-x f (t+dt)]P(f,t+dt)P(r,t\x,t), 

(6.8) 
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where P(r,t\x,t) is defined in Eq. (5.9) and Xf(t + dt) is 
given by Eq. (5.7). For the system we are considering, to 
find A(x,t) and D(x,t) we simply difïcrentiatc Eq. (6.1) 
and equate this with Eq. (5.7). Doing this gives 

A(x,t) = -Bx + Bm(t), (6.9) 
D(x,t) = B. (6.10) 

Combining the quantum and classical parts of the op- 
eration and using the same techniques as in Sec. VB 
allows us to rewrite Eq. (6.3) for continuous-in-time mea- 
surements as 



dp R (x,t) = dt(Bd x x + ±B 2 d£+£)p R (x,t) 



-dt 



(xn)tMt + dt) - (x R ) t 



-dt^d x B[&p R (x,t) + pa.(x,t)&i], (6.11) 
where (x R )t = J xTr[p R (x, t)]dx and 

r(t + dt)dt= (x R ) t dt+ y/f3dW{t). (6.12) 

This cquation (6.11) has been labeled the Superoperator- 
Kushner-Stratonovich equation [22] and represents the 
evolution of the combined supersystem. The first line 
contains the free evolution for both the quantum and the 
classical systems. For this quantum system 



C[à]p 



-in 



a x ,p] +jV[&}p, 



(6.13) 



where Í7 is the Rabi frcqucncy and T> is the damping su- 
pcropcrator and is defined in Eq. (2.14). The sccond line 
of Eq. (6.11) describes the gaining of knowledge about 
the state of classical system via Gaussian measurcmcnts. 
Lastly the third line describes the coupling of the quan- 
tum and classical system. 

For a TLA we can write the state of the supersystem 

as 

p(x,t) = í[P{x,t)Í + X(x,t)a x + Y(x,t)& y + Z(x,t)a z ]. 

(6.14) 

Note that P(x, t) is the marginal state of knowledge for 
the classical system (found via tracing out the quantum 
degrees of freedom). Substituting this into Eq. (6.11) 
gives the following four couplcd partial diffcrcntial cqua- 
tions 

Pr = [x- (aa) t ][r- (x R ) t }P R /(3 + Bd x [xP R 

-V7*r] + \B 2 dlP R (6.15) 
X R = [x - (x R ) t ][r - (x R ) t ]X R /P + \B 2 dlX R 

+Bd x [xX R - ^P R - ^fZ R ] - í 7 Ar (6.16) 

Y R = [x-(x R )t}[r~(x R ) t }Y R /f3 + Bd x [xY R } 

+\B 2 dlY R - nz R - i 7 r R (6.17) 

Z R = [x-(x R )t)[r-(x R ) t ]Z R /(3 + Bd x [xZ R 

+V7 X R ] + \B 2 d 2 x Z R + QY R - 1 (P R + Z R ). 

(6.18) 




FIG. 8: pn(x,t) calculated via numerical integration (sòlid) 
and via the ostensible method for an ensemble size of 10 000 
(dotted). The parameters are (3 = 0.5, B = 2, 7 = 1 and 
Q. — 5 and initial conditions p(x,0) = P{x)\g)(g\ where P(x) 
is a Gaussian with mean zero and variance 0.1. 



To determine the state of knowledge for the quantum 
system we simply integrate out the classical degrees of 
freedom. 

To illustrate a trajectory for this supersystem the fol- 
lowing parameters were used; /3 = 0.5,S = 2,7=1 
and Q = 5. The results are shown in Fig. 8 (sòlid line) 
for a randomly chosen record R. This figurc displays the 
mean and the variance of the classical trajectory found 
via tracing over the quantum degrees of freedom as well 
as the quantum state in Bloch representation after we 
have integrated out the classical degrees of freedom. 



B. Fictitious trajectòries: The ostensible numerical 
technique 

In the above section we observed that to be able to 
calculate the supersystem trajectory we needed to solve 
four coupled partial differential functions (each involving 
derivatives with respect to a classical configuration coor- 
dinatc x). This is a rather lengthy calculation which for 
highcr dimcnsional (d) quantum systems will require d 2 
partial diffcrcntial equations. Here we present our linear 
method that allows us to reduce the problem to d + 2 
couple differential equations. The expense, again, is that 
an ensemble average must be performed. 

To do this we simply note that we can define the fol- 
lowing quantum and classical states 



Of(t + dt,t)p(t) 
p f (t + dt) = n J/£ \ (6.19) 



A(/) 



and 



p , / JdxÓ rJ (x',t + dt\x,t)P(x,t) 
P rJ (x,t + dt) = A(/)A(r) ' ( } 



13 



whcrc 

6 r f (x',t + dt\x,t) = 5[x' - x f (t + dt)]P(r,t\x,t)A(f). 

(6.21) 

Note the bar above A(/) means that the ostensible dis- 
tribution used to scale the quantum state does not have 
to be the same as that used to scale the classical state. 
Here for simplicity we consider only the case when they 
are the same (as no numerical advantage is gain by differ- 
ent choiccs). Using the above equations we can rewrite 
Eqs. (6.3) and (6.5) as 



p r {x' ,i + dt) = 



P(r, t + dt) 



E/ Mf)Mr)Prj(x',t + dt)p f (t + dt) 
P{r, t + dt) 



dx 



f 



(6.22) 

Tr[A(f)A(r)P r ,f{x,t + dt) 



xp f (t + dt)}. 



(6.23) 



Thus to simulate pn(x,t) we nccd only to calculate 
Pn,F(x,t) and p~F(t) for a specific record R. 

For the above TLA-classical detector system with A(r) 
and A(/) defined by Eqs. (5.15) and (4.15) respectively, 
Pr,f{x', t) has a solution of the form pR,(t)S[x' — £r,f(í)] 
where xf(í) is given by 



dx F (t) = dt[-Bx F (t) + Bf(t + dt)} 
and p-R.,F(t) is given by 

<*r,f(í) = dt[x F {t) - A][r(í + dt) - A]p R)P (í)//3, 
Thus we can rewrite Eq. (6.22) as 

Pn(x,t) 



(6.24) 



(6.25) 



(6.26) 



where p F (t) = Ti[p F (t)]. 

To determine the evolution of the ostensible quantum 
state we simply substitute the measurement operator de- 
fined in Eq. (6.6) with H = ilà x /2 and the ostensible 
distribution A(/) into Eq. (6.19). Doing this gives 



Ep 


S[x- 


£f(í)W,f(í)pf(í) 


E F 


PR,F(,t)p F (t) 





dp F {t) 



-iQ r „ . . . 
dt— —[<J x ,PfW\ 



2 

dt[f(t + dt) 
+V7PF(í)ò· t - pp F {t)] 



dt-fV[à]p F (t) 
P\[\/ïèpF{t) 



(6.27) 



Howcver since we have assumed that all the quantum sig- 
nal is fed into the detector the evolution of the ostensible 
quantum state can be written as an ostensible SSE. That 
is, 

d\Mt)) = dt(- t —à x + [f(t + dt)- l· í}(^jà- l· L/2) 
-hh^a - ^jfià + p 2 /4}) \Mt)). (6.28) 



Thus to determine pr(x,í) all we nced to do is solve 
the above SSE and Eqs. (6.24) and (6.25) for R assumed 
known and f(t + dt)dt = dW + dtp, where dW is a Wiener 
increment. Once solved the quantum state conditioned 
on R is given by 



Xn(t) 



Et 



E t 



Pr,f(í)Xf(í) 



í»r,f(í)pf(í) 



(6.29) 



where \i = { x i, Vi, ii\ are the Bloch vectors of the quan- 
tum state. The moments of the classical state are given 
by 



Ep 


Xp(í)pr,f(í)pf 


(*)" 


E F 


Pr,f(í)íf(í) 





(6.30) 



To illustrate this method we considered two choices 
for the ostensible distributions. The first is A = fi = 0; 
that is, all the ostensible distributions arc Gaussian 
distributions of mean zero and variance dt. The sec- 
ond case corresponds to the situation when A = and 
/' = y/l^Hà + & )Pf(*)]í that is, the fictitious distri- 
bution is treated as the real unobservable distribution. 
Both cases were simulated to show the robustness of our 
numerical tcchnique and to demonstrate that while any 
ostensible distributions can be chosen a more realistic 
choice will result in a faster convergence. To demon- 
strate this we numerically solved Eq. (6.11) and used 
this as our reference solution. Then we compared the 
mean and variance of the classical marginal states and 
the fidelity for quantum reduced states (using Eq. (4.21) 
once the classical space has been removed) for both os- 
tensible cases and with ensemblc sizes of 100 and 10 000. 
These results are shown Fig. 9 where it is observed that 
for the larger ensemble size the difference in the classical 
marginal state is small and the quantum fidelity is close 
to one, indicating that our ostensible method has repro- 
duced the known result and is converging. Furthcrmorc 
it is observed that for the second case for the same en- 
semble size this difference is smaller thereby indicating 
that the second method convergence is faster. 



VII. DISCUSSION AND CONCLUSION 

The central tòpic of this paper was to investigate the 
conditional dynamics of partially observed systems (clas- 
sical and quantum) . Due to the fact that the information 
obtained is incomplete we have to assign a mixed state to 
the system. For a quantum system this means the state 
of knowledge given result r is given by the state matrix 
p r (t) and for a classical system a probability distribution 
P r (x, t) has to be used. If we consider a joint system (for 
example a classical detector is used to monitor a quantum 
system) the conditional state is given by p r {x,t). 

Even when we consider continuous-in-time monitoring 
we can still have incomplete information because of un- 
observed processes. For this case the conditional state 



14 



A 2 
X 
. 2 


A 






















„ 2 

~S 
> 

to 2 


A 


- . t^>-^ — t ï » ^^tjw- 




















C 1 

O 

0.9 


A 
































K 1 

1 : 


B 




















to - 1 






„ 2 

~í 
> 

to 2 


B 






















~ 1 

3 

u 


B 










"■"=0 1 2 3 4 £ 



time, t [y 1 ] 

FIG. 9: This figure shows the quantum and classical fidelity 
between the actual solution and our ostensible solution for 
pn.(x,t). Part A corresponds to the X — p = case while 
part B represents the A = and p — ^/7Tr[(<7 + à^pF^)] 
case. In both case an ensemble size of n — 100 (dotted) and 
n = 10000 (sòlid) was used. The system parameters are the 
same as in Fig. 8. 



trajectòries obey either a stochastic master equation (for 
a quantum system), a Kushncr Stratonovich equation 
(for a classical system) or a superoperator Kushncr- 
Stratonovich equation (for the joint system). That is, to 
simulate the conditional state we have to solve a rather 
numcrically cxpcnsive equation. In this paper we showed 
that by introducing a fictitious record F for the unob- 
served processes and ostensible measurement theory we 
can reduce this problem to solving pure states (stochas- 
tic Schròdinger equations for the quantum system or 
stochastic diffcrcntial equations for the classical system) 
conditioned on both R and F. Then by averaging over 
all possible F we get the require conditional state. That 
is the numerical memory requirements are decreascd by 
a factor of N, the number of basis states for the system. 
However, this is at the cost of an ensemble average. 

In summary, our ostensible method will be uscful for 
investigating realistic situations where the dimensions of 
the systems are large. It is also much easier to implement 
numcrically than the Standard technique, so we expect it 
to find immediate applications. 
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APPENDIX A: WHY IT IS NECESSARY TO USE 
THE OSTENSIBLE METHOD 

To show that we must usc an ostensible distribution 
rather then the real distribution, for our numerical tech- 
nique, we consider two consecutive measurements. The 
state of the system (which we take to be quantum for 
specificity) after these two measurement is 

Pr 2 ,rx = ^2 P(f2,fl\r2,r 1 )pr 2 J 2 ,r 1 J 1 



P{ri,r 2 ) 



(Al) 



This can be rewritten as 



PT2 ' T1 fi P(r 2 ,/2|r 1 ,/ 1 )P(r 1 ,/ 1 ) P(r u r 2 ) ' 

J2-J1 

(A2) 

The first term can be viewed as the part that determines 
the trajectory and the second term as the part which 
determines the weighting factor for this trajectory. Con- 
sidering only the weighting factor we can rewrite this as 
P(Í2, fi 1^2, t'í), which, unless we have the full numerical 
solution, is not determinable. To be more specific we can- 
not separate this term into P (fi)P(f 2 \fi, fi) and thus we 
cannot create a trajectory that steps through time with 
the correct statistics for fk- 

However by introducing an ostensible distribution we 
can rewrite Eq. (A2) as 



Pr 2 ,ri 



Pf2,r 2 Jl,ri~ 



A( / 2 ,r 2 ,/i,ri) 
P(ri,r 2 ) : 



where 



Pf2,r 2 Ji,r 1 



J[Mf 2 ,r 2 M fl ,rM0) 

Hf2,r 2 \fi,n)A(f 1 ,n) 



(A3) 



(A4) 



and we have completo freedom to choose any 
A(/2, r 2 , /i, ri). As such we are not restricted to using 
the undcterminablc distribution P(f 2 , r 2 , /i, r\). This 
implics that to unravel the conditional state conditioned 
on some real record R in terms of hetitious results the 
corresponding trajectòries must be unnormalized. 

In this paper we make two choices for 
A(/ 2 ,f 2 ,/i,ri). The first choice is A(/ 2 , r 2 , /i, n) = 
A(/ 2 )A(/ 1 )A(r 2 )A(n) where A(r k ) and A(/ fc ) are 
Gaussian distributions of variance dt and mean 
A and p respectively. The second choice is 
A(/ 2 ,r 2 ,/!,n) = A(/ 2 |/i,n)A(/i)A(r 2 )A(n) where 
A(rk) is the same as beforc but the fictitious results 
were chosen based on the true probability we would 
expect based on the past real and fictitious results 
up to but not including the current time. That is 
P = V^Tr[x 2 p R , F (t)]/Tr[p RiF (í)]. 

A third possible choice would be A(/ 2 , r 2 , fi, ri) = 
A(/ 2 |/i,ri)A(/i)A(r 2 |ri,/i)A(ri), where both the real 
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FIG. 10: This figure shows a typical trajectory for a radia- 
tively damped TLA undergoing homodyne-a; detection of effi- 
ciency 0.4 (sòlid line). The dotted line represents the solution 
using our method whereas the dashed line represents the BG 
method. Both our method and the BG method were imple- 
mented using an ensemble size of 10000. 
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FIG. 11: The top two plots show the difference between the 
exact solution and the BG method for n = 1000 and 10000. 
The bottom two plots show the difference between the exact 
solution and our method for n = 1000 and 10000. 



and fictitious distribution are chosen based on the 
past results up to but not including the current time. 
That is A = y /^nTï{± 1 p R ^ F {t)\/i:r[p R . F {t)] and A* = 
^/T^Tr[x2PR, i p(t)]/Tr[pR i p(í)]. This would seem the clos- 
est choice to P{f2, ?"2, fi, fi) but it is important to note 
that it is not the same. It is still an ostensible distri- 
bution, because the true distribution is based upon the 
entire measurement record, including results in the fu- 
ture. Thus our trajectory equations will still bc unnor- 
malized and we cannot replace (r — X)dt with dW{t) as 
Tr[fipR iP (í)]/Tr[pa,F(í)] is not equal to Tr[xip R (t)] for 
all possible fictitious records. If we were to make this sub- 
stitution we would generate normalized equations, but 
averaging over all possible fictitious record would not give 



a typical trajectory for the state conditioned on the par- 
tial record R. This is precisely the mistake Brun and 
Goan (BG) make in reference [20]. 

To be more specific let us consider their approach 
and our approach for the following simple system: A 
two level atom radiatively damped and monitored using 
homodyne-cc detection with an efficiency 77. That is this 
system is described by the SME 

dp R (t) = dW[à]p n + ^}H[a}p R dW(t), (A5) 

where r(t + dt)dt = dW(t) + dt^/fj (à x ). Now using BG's 
theory we would extend this equation to 

dp RF {t) = dW[à]p RlF + y/rjH[à]p RtF dW(t) 

+ ^T^7,H[à]PK.FdW(t), (A6) 

which has a pure state solution. Now they argue that 
by ensemble averaging over the fictitious noise process 
dW Eq. (A5) is recovered. However this is incorrect. 
Although 

v w [Vvn[à}p R . F dw(t)} = Vv^w[n[à]p R . F }dw(t), 

(A7) 

the nonlinearity in the superoperator TL means that 

E w [H[à}p RF ] ï H[à}E w [pR,F] = H[à}p R . (A8) 

To see this explicitly we expand Eq. (A6) for two mea- 
surement (two steps in time) 

Pri,r2,fi,f 2 ( 2 dt) = 

po + dtv[d] PriJl + ^nH[à]p ruh dW 2 {t) + 

^1 - T 1 H[à}p ri ,f 1 dW 2 (t) + dtt>[à]p + 
y/rjH[8]podWi(t) + yï^HWpodW^t). (A9) 
Looking at this equation we see that the problem term is 

^ h [H[à]p riJl \ = E fl [àp riJl + p riJl à 
-Tr[à x p rufl ]p rufl ] 
^ àp ri + p ri a - Tr[à x p ri ]p ri (A10) 

due to the non-linearity. 

In our method this problem does not oceur because 
we use ostensible distributions and linear equations. For 
this system the unnormalized state is 

dp R , F (t) = dW[à]p R . F + y/rjÜ\[à]p RtF dt(r - X) + 
y/Ï^Hpl&lpRjrdWtt), (All) 

where 7í x [À}p = Ap + pA — \p with \ being either A or 
p. If we expand this to two measurements we find that 
the problem term does not oceur as 

E /i [#[o-]Pri ,/i ] = E/, [Òp ri ifl + p ri ifl <T - Xp ri ,/, ] 

= &p ri + p ri ò - Xp ri . (A12) 

To show the magnitude of the error in BG's ap- 
proach we numerically solved the SME for 77 = 0.4 us- 
ing Eq. (A5), BG's method and our method. We first 
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randomly gcncrate and store a string of dW's, and use 
these to generate the true record R via Eq. (A5). For 
BG's method we use this string of dW's, and generate an 
ensemble using randomly generate strings of dWs. The 
resultant solution would, according to BG, correspond to 
the solution of Eq. (A5). For our method, wc use the true 
record R and again randomly generate strings of dWs to 
obtain an ensemble average. The rcsults of these simula- 
tions are shown in Figs. 10 and 11. Hcrc wc sec that BG 



thcory disagrees significantly with the exact rcsult. Fur- 
thermore this discrepaney cannot be a statistical error 
as Fig. 11 shows that when the ensemble average is in- 
creased this difference remains approximately constant. 
That is, this simulation confirms that BG's method fails 
to reproduce a typical solution to Eq. (A5). By contrast 
our method reproduces Eq. (A5) to within statistical er- 
ror. 
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